Quantum correlations in spin chains at finite temperatures and quantum phase 

transitions 
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We compute the quantum correlation (quantum discord (QD)) and the entanglement (EoF) be- 
tween nearest neighbor qubits (spin-1/2) in an infinite chain described by the Heisenberg model 
(XXZ Hamiltonian) at finite temperatures. The chain is in the thermodynamic limit and thermal- 
ized with a reservoir at temperature T (canonical ensemble). We show that QD, in contrast to EoF 
and other thermodynamic quantities, spotlight the critical points associated to quantum phase tran- 
sitions (QPT) for this model even at finite T. This remarkable property of QD may have important 
implications for experimental characterization of QPTs when one is unable to reach temperatures 
below which a QPT can be seen. 
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Quantum phase transition (QPT) is a purely quan- 
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turn process [l[ occurring at absolute zero temperature 
(T = 0), where no thermal fluctuations exist and hence 
no classical phase transition is allowed to occur. QPT 
is caused by changing the system's Hamiltonian, such 
as an external magnetic field or the coupling constant. 
These quantities are generally known as the tuning pa- 
rameter. As one changes the Hamiltonian one may reach 
a special point (critical point) where the ground state 
of the system suffers an abrupt change mapped to a 
macroscopic change in the system's properties. This 
change of phase is solely due to quantum fluctuations, 
which exist at T = due to the Heisenberg uncer- 
tainty principle. This whole process is called QPT. 
The paramagnetic-ferromagnetic transition in some met- 
als 0, the superconductor- insulator transition Q, and 
superfluid-Mott insulator transition are remarkable 
examples of this sort of phase transition. 

In principle QPTs occur at T = 0, which is unattain- 
able experimentally due to the third law of thermody- 
namics. Hence, one must work at very small T, as close 
as possible to the absolute zero, in order to detect a QPT. 
More precisely, one needs to work at regimes in which 
thermal fluctuations are insufficient to drive the system 
from its ground to excited states. In this scenario quan- 
tum fluctuations dominate and one is able to measure a 
QPT. 

So far the theoretical tools available to determine the 
critical points (CP) for a given Hamiltonian assume 
T — 0. For spin chains, for instance, the CPs are de- 
termined studying, as one varies the tuning parameter, 
the behavior of either its magnetization, or bipartite Q 
and multipartite @ entanglement, or its quantum cor- 
relation (QC) 0. By investigating the extremal values 
of these quantities as well as the behavior of their first 
and second order derivatives one is able to spotlight the 
CP. However, the T = assumption limits a direct con- 
nection between these theoretical "CP detectors" and ex- 



periment. Indeed, if thermal fluctuations are not small 
enough excited states become relevant and the tools de- 
veloped so far cannot be employed to clearly indicate the 
CP. 

In this Letter we remove this limitation and present 
a theoretical tool that is able to clearly detect CPs for 
QPTs at finite T. We show that the behavior of strictly 
QCs at finite T, as given by the thermal quantum 
discord (TQD) 0, unambiguously detects CPs for QPTs 
that could only be seen, using previous methods, at T = 
0. This remarkable property of TQD, on one hand, is 
an important tool that can be readily applied to reduce 
the experimental demands to determine CPs for QPTs; 
or even allow such a detection for those systems where to- 
day's technology makes it virtually impossible to achieve 
the necessary T below which quantum fiuctuations dom- 
inate. One the other hand, this characteristic of TQD 
shows that QPTs have a decisive influence on a system's 
physical property not only for small T but also above T 
where quantum fluctuations no longer dominate. 

In order to show that TQD detects a QPT at finite 
T, we study the anisotropic spin-1/2 Heisenberg chain 
(XXZ) in the thermodynamic limit. We assume the infi- 
nite chain to be in thermal equilibrium with a reservoir 
at temperature T, i.e., its density matrix is described 
by the canonical ensemble. Tracing out all spins but 
the two nearest-neighbors we get their reduced density 
matrix as a function of two-point correlation functions, 
which are evaluated by sol ving a set of non-linear inte- 
gral equations (NLIE) [lol . Ill| . The two-qubit density 
matrix allows us to compute TQD and investigate its 
properties for T > as we change the system's Hamilto- 
nian. We show that TQD is maximal and its first order 
derivative with respect to the tuning parameter is dis- 
continuous at the quantum CP, not only at T = 0, 
but also at T > 0. This behavior is robust enough to be 
seen for high T. Furthermore, we have also computed the 
entropy, magnetization, magnetic susceptibility, and spe- 



2 



cific heat, for the whole chain, and two-site correlations 
between the two nearest-neighbor spins as well as their 
entanglement. We show that none of these quantities de- 
tect unambiguously the CP for T > 0. We also discuss 
why TQD possesses such a unique behavior in contrast 
to another quantity, namely, the entanglement between 
the two nearest-neighbors. 

The XXZ Hamiltonian can be written as 



(1) 



Note that EoF and QD coincide for bipartite pure states; 
for mixed states, though, their difference becomes mani- 
fest being both zero, however, when only classical corre- 
lations are present. We can also conceptually understand 
QCs in comparison with entanglement by noting that the 
latter is due to the superposition principle applied to the 
whole Hilbert space of a bipartite system. However, QCs 
as given by QD captures, on top of that, the correlations 
coming from superposition of states within each subsys- 
tem, a purely quantum effect that it is not possible clas- 
sically [16[. From this perspective, one can better grasp 
why there exist states with zero entanglement but finite 



where periodic boundary conditions are assumed and A QCs [17|. Another interesting and operational interpre- 
is the anisotropy parameter. Here L — > 00 and crjjcrj, 
and (t| are the usual Pauli matrices acting on the j-th 
qubit. Throughout this Letter H = I and J = 1 unless 
noted otherwise. At T = the XXZ model has two CPs 
At A = 1 we have a continuous phase transition 
and at A = — 1 we have a first-order transition. 

The density matrix for a system in equilibrium with a 
thermal reservoir is p ~ exp {~j3H)/Z , where j3 = 1/kT, 
Z = Tt {exp {—(3H)} is the partition function, and the 
Boltzmann's constant k is set to unity. The nearest- 
neighbor two qubit state is obtained by tracing all but 
the first two spins, pi2 = Tri_2{p}. Due to the transla- 
tion invariance and U{1) invariance {[H,J2j=i<^j] = 0) 
of ID), we can write the reduced density matrix as follows. 



tation for QD is achieved looking at the thermodynamic 
properties of a quantum system. In [isj it is shown that 
QD is related to the difference of work that can be ex- 
tracted acting either globally or locally at a heat bath 
with a bipartite state when one-way communication is 
allowed. 

For state ^ QB is ^ QD = [g{l - 2dx - d^) 
+ 2 ,9(1 + d,) + .9(1 + 2d, - 4)] /4 - [.g(l + D) + 
g{l - D)] /2, with d^ = (cr^crf ), d, (crfcr?), and 



Z? = Max{|Kaf)|,|Ka|)|}, 



(5) 
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(2) 



These two-point correlation functions can be written in 
its simplest form in terms of derivatives of the free energy 
/ = (-l//3)limL^oo(lnZ)/L, 

(cT|a|+i) = dAf/J, (^J^J+i) = (" - A9a/)/2J,(3) 

with u = dp{Pf) the internal energy. In order to deter- 
mine the free energy in the thermodynamic limit and at 
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finite T one has to solve a suitable set of NLIE 

Now we can use ^ in order to show that the entan- 
glement, as measured by the entanglement of formation 
(EoF) Q, is^oF_= -g{f{C)) - g{l - f{C)), with 
/(C) = (1 + Vl~C^)/2, g{f) = /log2(/), and 



C Max{0, \{afa^ 



K^2)l/2} 



(4) 



the concurrence, an entanglement monotone. EoF quan- 
tifies a class of QCs that cannot be created by local op- 
erations and classical communication (LOCC) only jlSf . 
Recently, however, it became clear that there exist more 
general QCs if one removes the LOCC restriction. These 
correlations are measured by the quantum discord (QD) 
and it is believed that QD quantifies all correlations 
between two systems that has a pure quantum origin. 



Note that either | (erf erf ) | or | (erf erf) | is responsible for the 
value of D. As will be seen, it is the interplay between 
these two correlations that is relevant in our understand- 
ing of why QD detects a QPT at finite T and EoF does 
not 0. 

We are now in a position to present the behavior of 
TQD and EoF between two nearest-neighbor qubits in 
an infinite spin chain at finite T. We first plot TQD and 
EoF, for several T, as a function of the tuning parameter 
A. This allows us to prove the main claim in this Letter, 
namely, that TQD detects a CP of a QPT at finite T 
while EoF does not. Looking at Fig. [T]we see that EoF 
is maximal in the CP A = 1 only at T = 0, agreeing 
with the results of 0]. As we increase T the maximum 
no longer occurs at A = 1, moving to the region where 
A > 1. Also, the higher T the farther from the CP is lo- 
cated the maximum of EoF. On the other hand, TQD is 
maximal at A = 1 when T = and does not appreciably 
move away for T < 3. Moreover, its first order derivative 
is discontinuous at the CP not only at T = but also at 
T > 0, a remarkable result showing that TQD inherits 
at T > all of its important properties previously seen 
only at T = 0. This discontinuity of the first derivative 
of TQD at A = 1 is our CP detector for non null T. 
In order to prove this unique behavior of TQD, we have 
computed for several T many thermodynamic quantities 
for the infinite spin chain and also the pairwise correla- 
tions as a function of the tuning parameter A. As can be 
seen in Fig. [21 none of these quantities can clearly detect 
the CP at T > 0. 

Due to subtleties of the NLIE at A = -1 (J > 0), it 
is convenient to investigate how TQD behaves near the 
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FIG. 1. (Color online) EoF (top) and QD (bottom) as func- 
tions of the tuning parameter A for the XXZ model in the 
thermodynamic limit. The inset depicts QD for high T near 
the CP. T increases from top to bottom. The curves for T = 
and T = 0.01 cannot be distinguished from the T = 0.1. Here 
and in the following graphics all quantities are dimensionless. 



namely, being discontinuous at the CP Q. 
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FIG. 3. (Color online) EoF and QD for a chain of 8 and 
10 qubits described by the XXZ model. QD detects both 
quantum critical points at finite T while EoF does not. 




In order to complement our results, we fix the 
anisotropy parameter at A = 1 and then vary the cou- 
pling constant J from negative to positive values, i.e., we 
go from a ferromagnetic to an antiferromagnetic regime. 
As can be seen in Fig. |4] TQD decreases as one varies J 
towards zero from both sides Q . Similar to the previous 
case, TQD inherits for finite T its behavior at T = 0. 
However, EoF is only non zero for the antiferromagnetic 
regime; and for finite T this only occurs away from the 
vicinity of J = 0. In other words, the behavior of TQD 
around J = and T > are qualitatively similar to its 
behavior at T = while this is no longer true for the 
behavior of EoF. 



FIG. 2. (Color online) Thermodynamic quantities for the 
XXZ model in the thermodynamic limit. The T = and 
T = 0.1 curves for the two-point correlation functions are 
indistinguishable. Note that at T = the magnetic suscepti- 
bility also detects the phase transition being discontinuous at 
the CP [l^]. The specific heat and entropy are null at T = 0. 



CP A = — 1 by means of numerical diagonalization of the 
Hamiltonian ([T]) [i^. We computed its thermal density 
matrix, and then calculated the nearest-neighbor reduced 
density matrix for lattice sizes L = 8 and 10. Again, only 
TQD was able to detect the CP for T > 0. Looking at 
Fig. [3]we clearly see that TQD successfully picks the CPs 
at A = ±1 while EoF does not. For finite T, the first 
derivative of TQD is discontinuous at both CPs. EoF, on 
the other hand, is zero around A = — 1 and its maximum 
gets shifted to the right at A = 1. Note that for small T 
and A = — 1 TQD also resembles its behavior at T = 0, 
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FIG. 4. (Color online) EoF and QD for a chain of 8 and 10 
qubits described by the XXX model. 

We can understand this unique aspect of TQD, spe- 



4 



cially in contrast to EoF, by taking a careful look at 
the analytical expressions giving EoF and TQD. The 
main difference in behavior between EoF and TQD is 
connected to Eqs. ^ and (O, being directly related 
to the maximization process leading to these quantities. 
For the XXZ model and at finite T, one can show that 
around the two CPs the function maximizing (j4|) does not 
abruptly change. It is either or |(tTfCTf)| — |l + (CTfCT|)|/2. 
On the other hand, for ([SJ, the function maximizing it 
changes exactly at the CPs. Before the CPs one has ei- 
ther |(crJ^(T|)| or Kcr^crf)! as the maximum but after them 
this role is exchanged. Indeed, in the vicinity of A < —1 
D is given by |(crfcr|)| while for —1 < A < 1 it is deter- 
mined by |((Ticrf)| (see Fig. [2]). Finally, in the vicinity 
of A > 1 it is determined by | (erf erf) |. It is this change 
in the function maximizing D, which occurs at T = (7[ 
and shown here also to occur at T > 0, that is responsi- 
ble for the discontinuity of the first derivate of TQD. For 
the XXX model, |(crftTf)| = |(criCT|)|, and therefore no 
cusp-like behavior for TQD is observed. However, TQD 
is only zero at J = for any T while EoF is always zero in 
the vicinity of J = for T > 0. Moreover, working with 
small chains (up to 10 qubits) for various T, we observed 
that the second derivative of TQD possesses a relatively 
high value near J = 0. We believe that it is likely that 
as one approaches the thermodynamic limit the peak of 
the second derivative moves towards J = 0. 

In summary, we presented a remarkable characteristic 
of quantum correlations as given by the quantum dis- 
cord (QD): its ability to detect critical points (CP) of 
quantum phase transitions (QPT) at finite T. Indeed, 
by solving an infinite chain described by the XXZ model 
in the thermodynamic limit, we showed that QD is able 
to highlight the CPs of QPTs for T > while neither the 
entanglement nor any thermodynamic quantity achieve 
the same feat. This property of QD may be useful in 
the experimental detection of CPs for QPTs where one 
is not able to reach the temperatures below which a QPT 
can be seen. Conceptually, this capacity of QD to detect 
CPs of QPTs for T > and its interes ting and puzzling 
dynamical robustness against noise |2J, [25| illustrate the 
broad range of scenarios where QD helps in the under- 
standing of fundamental issues of quantum mechanics. 

TW thanks CNPq for funding. CT thanks Volkswa- 
gen Foundation for financial support. GAPR thanks 
FAPESP for funding. GR thanks CNPq/FAPESP for 
financial support through the National Institute of Sci- 
ence and Technology for Quantum Information. GAPR 
and CT thank Frank Gohmann for discussions. 



* |rigolin@ufscar.br] 

[1] S. Sachdev, Quantum Phase Transitions (Cambridge 
University Press, Cambridge, 1999). 



[2] S. Rowley et ai, Phys. Status Solid B 247, 469 (2010). 
[3] V. F. Gantmakher and V. T. Dolgopolov, Physics- 

Uspekhi 53, 1 (2010). 
[4] M. Greiner et ai, Nature (London) 415, 39 (2002). 
[5] L.-A. Wu, M. S. Sarandy, and D. A. Lidar, Phys. Rev. 

Lett. 93, 250404 (2004). 
[6] T. R. de Oliveira et ai, Phys. Rev. Lett. 97, 170401 

(2006); T. R. de Oliveira et al., Phys. Rev. A 73, 

010305(R) (2006); T. R. de Oliveira et al, Phys. Rev. 

A 77, 032325 (2008). 
[7] R. Dillenschneider, Phys. Rev. B 78, 224413 (2008); M. 

S. Sarandy, Phys. Rev. A 80, 022108 (2009). 
[8] H. OUivier and W. H. Zurek, Phys. Rev. Lett. 88, 017901 

(2001); B. P. Lanyon et al, ibid. 101, 200501 (2008); L. 

Henderson and V. Vedral, J. Phys. A 34, 6899 (2001). 
[9] T. Werlang and G. Rigolin, Phys. Rev. A 81, 044101 

(2010). Note that TQD refers to QD computed for states 

described by the canonical ensemble. 
[10] A. Kliimper, Z. Phys. B 91, 507 (1993). 
[11] M. Bortz and F. Gohmann, Eur. Phys. J. B 46, 399 

(2005); H.E. Boos et al, J. Stat. Phys., P08010 (2008); 

C. Trippe, F. Gohmann, and A. Kliimper, Eur. Phys. J. 

B 73, 253 (2010). 
[12] M. Takahashi, Thermodynamics of one-dimensional solv- 
able models (Cambridge University Press, Cambridge, 

1999). 

[13] Details of the calculations leading to the two-point cor- 
relation functions and a complete analysis of this model 
with and without magnetic field will be given elsewhere. 

[14] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998); T. Yu 
and J. H. Eberly, Quantum Inf. Comp. 7, 459 (2007). 

[15] R. F. Werner, Phys. Rev. A 40, 4277 (1989). 

[16] The superposition principle allows a subsystem of a bi- 
partite system to be prepared in a superposition of two 
non-orthogonal states. In such a case the whole system 
cannot be described classically resulting in QD 0. 

[17] K. Modi et al, Phys. Rev. Lett. 104, 080501 (2010). 

[18] M. Horodecki, J. Oppenheim, and R. Horodeckil, Phys. 
Rev. Lett. 89, 240403 (2002); I. Devetak, Phys. Rev. 
A 71, 062303 (2005); M. Horodecki et al., Phys. Rev. 
A 71, 062307 (2005); W. H. Zurek, Phys. Rev. A 67, 
012320 (2003); R. Dillenschneider, Eurphys. Lett. 88, 
50003 (2009). 

[19] S. Luo, Phys. Rev. A 77, 042303 (2008). 

[20] Technically, which correlation function dominates is re- 
lated to the maximization, over all possible measure- 
ments on one of the subsystems, carried out to compute 
an extension of the classical mutual information (MI) for 
quantum states using the Baye's rule. This quantity is 
then subtracted from another way of defining MI without 
Baye's rule. The difference between these two quantities 
leads to QD [^. For states as ((2]) the maximization lead- 
ing to QD is achieved by either one of two different sets 
of measurements, defining which correlation function is 
the relevant one in the expression for QD [2ll]. 

[21] J. Maziero et al., Phys. Rev. A 80, 044102 (2009); F.F. 
Fanchini et al., Phys. Rev. A 81, 052107 (2010). 

[22] C.N. Yang and CP. Yang, Phys. Rev. 151, 258 (1966). 

[23] It is worth noting that the two-point correlations for L = 
10 are very close to those (Fig. [2)) at the thermodynamic 
limit for A > 0. 

[24] L. Mazzola, J. Piilo, and S. Maniscalco Phys. Rev. Lett. 
104, 200401 (2010). 

[25] Jin-Shi Xu et ai, Nat. Commun. 1 (2010) 7. 



